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Abstract 

We present a method to characterize the polarization state of a hght field in the 
continuous-variable regime. Instead of using the abstract formalism of SU(2) qua- 
sidistributions, we model polarization in the classical spirit by superposing two 
harmonic oscillators of the same angular frequency along two orthogonal axes. By 
describing each oscillator by a s-parametrized quasidistribution, we derive in a con- 
sistent way the final function for the polarization. We compare with previous ap- 
proaches and discuss how this formalism works in some relevant examples. 
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1 Introduction 

Polarization is a fundamental property of light, both in the quantum and in the 
classical domain. In quantum optics, polarization has been mainly examined in 
the single-photon regime [1,2,3,4,5,6,7]. Nevertheless, different schemes have 
been proposed [8,9] and experimentally implemented [10,11] to characterize 
the continuous-variable limit of the quantum Stokes parameters. We stress 
that these continuous-variable polarization states can be carried by a bright 
laser beam providing high bandwidth capabilities and therefore faster signal 
transfer rates than single photon systems. In addition, they retain the single- 
photon advantage of not requiring the universal local oscillator necessary for 
other proposed continuous variable quantum networks. 

Since the seminal paper of Wigner [12], and the major contributions of Moyal [13], 
Stratonovich [14], and Berezin [15], it seems indisputable that phase-space 
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methods, based on using quasidistributions that reflect the noncommutatibil- 
ity of quantum observables, constitute a valuable tool in examining continuous 
variables in quantum optics [16,17,18]. 

In particular, these methods have had great success in analyzing one-mode 
fields; i. e., Hciscnbcrg-Weyl quasidistributions representing the quantum dy- 
namics in the flat q-p (or, equivalently, a-a*) space. Although not so popular in 
the quantum-optics community, spinlike systems, with the sphere ^2 as phase 
space, have been also discussed at length in this framework [19,20,21,22,23,24,25,26,27]. 
The resulting functions, naturally related to the SU(2) dynamical group, have 
been used to visualize, e. g., nonclassical properties of a collection of two- level 
atoms [28]. 

Since the Stokes operators can be formally identified with an angular mo- 
mentum [29,30,31,32], one may naively expect a direct translation of these 
SU(2) quasidistributions to the problem of polarization. However, this is not 
the case, mainly because they act on different types of Hilbert space [33]. We 
can then conclude that the problem of an adequate quasiclassical description 
of polarization of light is still an open question [34,35]. 

The Stokes operators arc a particular case of the Schwinger map [36] , with two 
kinematically independent oscillators. In this spirit, it has been recently shown 
that the Stokes operators are the constants of motion of the two-dimensional 
isotropic harmonic oscillator [37]. This refiects the fact that the polarization 
of a classical field can be adequately viewed as the Lissajous figure traced 
out by the end of the electric vector of a monochromatic field [38]. In sharp 
contrast, in quantum optics the probability distribution for the electric field 
can be very far from having an elliptical form [39,40]. 

We wish to investigate this point from the perspective of quasidistributions. 
We show that one can start from the s-ordered quasidistributions from two 
kinematically independent oscillators: by eliminating an unessential common 
phase, we get well-behaved quasidistributions on the Poincare sphere. We 
apply the resulting family of polarization quasidistributions to some relevant 
states, and conclude that they constitute an appropriate tool to deal with such 
a basic variable. 



2 Phase-space representation of a hcirmonic oscillator 

To keep the discussion as self-contained as possible, we first briefly summarize 
the essential ingredients of phase-space functions for a harmonic oscillator that 
we shall need for our purposes. 
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In the Hilbcrt space 7i, the state of the system is fully represented by its 
density operator q. In the phase-space formalism, g is mapped by a family 
of functions (quasidistributions) W^^\a) onto the classical phase space X 
{a & X). This map is usually implemented by the generalized Weyl rule [26] 

W^'\a)^T^[QW^'\a)l (1) 
where the generating kernel w^^\a) fulfill the general properties 

(2) 

X 



The index s that labels functions in the family is related to the s ordering. The 
values —1, +1, and correspond to the normal, antinormal, and symmetric 
ordering, respectively, or equivalently to the P, Q, and W functions. We stress 
that these quasidistributions can be determined in practice by using simple and 
efficient experimental procedures [41,42,43,44,45,46]. Moreover, they provide 
a simple measure of the nonclassical behavior of quantum states [47,48,49,50]. 

Let us now turn to the outstanding case of a harmonic oscillator represented 
by annihilation and creation operators d and d^, that obey the canonical com- 
mutation relation 

[d,d^] = i. (3) 



The phase space is the complex plane C and the invariant measure is dn{a) ~ 
d'^a/n'^. The operator 

D{cy) — exp(Q;a^ — a*d) (4) 



is the standard displacement operator in the complex plane a and leads to 
introduce the standard coherent states as 

\a)^D{a)\0), (5) 



where |0) denotes the ground state. In this case, the kernel w^^\a) is exactly 
the Cahill-Glauber kernel 

= i- / d2/5e-^"^''/2 exp(a/3* - a*p) Dia). (6) 
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It is often useful to represent this operator in a somewhat different form. After 
some calculations, it turns out that (6) may be rewritten as [50] 



a) = Y^^D{a) [-^j DKal (7) 



which can also be represented in a disentangled form 

= ^ e ^-—J exp {-^/) cxp (^aj . (8) 

The reader is referred to e.g. Ref. [17] to sec how this s-parametrized repre- 
sentation works for some elementary field states. 



3 Phase-space description of two orthogonal harmonic oscillators 



In classical optics, the superposition of two oscillations of the same angular 
frequency ou, one along the horizontal axis H and the other along the vertical 
axis V, results in an elliptical motion. This is the simplest Lissajous figure 
and is the basic physics behind the notion of light polarization. 

To translate this picture into quantum optics, we assume that the two oscil- 
lators are represented by the complex amplitude operators an {alj) and oy 
{a\r), fulfilling (3); i. e., 

[aj,al]^Sjki, j,k^H,V. (9) 



In phase space, these oscillators can be appropriately described by the product 
of the corresponding kernel operators 

w^'\aH, av) = w^"\aH) w^"\av). (10) 



A Lissajous figure needs only three independent quantities to be fully charac- 
terized: the amplitudes of each oscillator and the relative phase between them. 
We therefore introduce the parametrization 

an ^re'^cos{e/2), ay = re'^e-''^ sm{e/2), (11) 



where 



= + layp (12) 
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is a radial variable related with the global intensity. The parameters d and 
can be interpreted as the polar and azimuthal angles, respectively, on the 
Poincare sphere 82- In terms of this parametrization, equation (10) can be 
recast as 



where 

N — OfjCiH + ayCiv 



2 ^2 + 



s- 1 



2r 
1-s 



(13) 
(14) 



is the operator representing the total number of excitations and 

U{9,<j))=exp ^- 



(15) 



We next proceed to integrate over the physically irrelevant global phase ( : 



N 



i=0 



£!2 



„2£ 



tlHOA). (16) 



Finally, we integrate over the radial variable r (with the weight 2r'^) to obtain 
the phase-space kernel over the sphere ^2 



^ Note, in passing, that this equation can be expressed compactly as 



1 - s ; V s - 1 



1 - s 



where : : means normal ordering and Jq denotes the Bessel function of first kind 
and zero order. 
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s + 1 
s-1 



N 



xu{e,<p) 



u\e,<p). (17) 



If we observe that 



^=0 



(18) 



we easily obtain 



^I;W(^,0) = ^ t/(^,0) 



s-1 



s - 1 



s + 1 

t A 



1-s 



1-s 



u\e,<t>). (19) 



This is our central and compact result, that we shall work out in detail in the 
rest of the paper. 

For the antisymmetric ordering (s = —1), which corresponds to the Q-function, 
(19) reduces to 



_|_ ]_ \ S|rav 
— -yj {d^H^'H + 1) 



u\eA). (20) 



To proceed further we notice that 



lim 



s + 1\"^ 



s-^-i Vs - 1 



N=0 



(21) 



is precisely the projector on the subspace with zero excitations in the vertical 
mode V. Henceforth, we shall denote by 



\N,k) := \N-k)H®\k)v 



(22) 



a state with N — k excitations in the horizontal mode H and k in the vertical 
mode V. In view of (21) we have 



w^-'\9, cf>)=J2iN + 1)\N, 9, 0)(7V, 9, cj>\ , 

N=0 



(23) 
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where \N,9,(j)) are the standard SU(2) coherent states [51] 



|7V,^» = C/(e,0)|7V,O) 

1/2 

/i\ N-k 




sin-j [cos-j e-'^^\N,k). (24) 
The Q function generated by this kernel reads as 

oo 

Q{e, <t>) = iv[^(-^)(^, 0)] = ^ (AT + 1) g(7v, <t>) , (25) 

N=0 

where 

Q{N,0,<f>)^{N,e,<l>\g\N,0,<f>). (26) 
The normahzation condition 

1 rfi]g(^,0) = i, (27) 

where dfl — sin OdOdcj) is the differential of solid angle, is automatically fulfilled. 

Equation (25) is a well-known result, which may be derived from a variety of 
methods. The point we wish to stress is here that (25) involves only diagonal 
elements between states with the same number of excitations. Because of the 
lack of the off-diagonal contributions of the form (A^, 9^ (j)\Q\N' , 9, 6) with N ^ 
N', the Q function takes the form of an average over the subspaces with 
definite total number of excitations. The role of the sum in N is to remove 
the total intensity from the description of the state. 

We next pass to the symmetrical order (s = 0), which corresponds to the 
Wigner function. Prom (19) we immediately get 

w^''\9,<P)^U{9,<f>) [(-l)»W(24a,j + l)l U\9,<P). (28) 



To properly transform the operator in brackets, it proves convenient to intro- 
duce the following operators 



Jy = ^(^hO'V - OyOi?) , (29) 
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which constitute the Schwinger map for two independent oscillators [36] and 
are the generators of the group SU(2). Since U{6,(f)) is an element of SU(2), 
it can be represented as [51] 

U{e, (p) = e"''^'^' e-'^^^ e''^•^^ (30) 



Next, we note that o^jjCiH = N/2 — and that 

u{e,<j>) lu^ie,<i))^j-n, (31) 



where J = {J^, Jy, Jz) and n = (sin 9 cos 0, sin 9 sin 0, cos 9) is a unit vector on 
the sphere «S2. Then, equation (28) simphfies to 

w^^\9, (j)) = {-lf^{N - 2J • n + i) U{9, 0) exp(-i7rJ^) U\9, 0). (32) 



This Wigner function will be examined for a variety of states in the next 
section. However, it seems pertinent to compare before these results with pre- 
vious approaches using the machinery of SU(2) quasidistributions. One can 
immediately check that both give the same Q function, meanwhile the corre- 
sponding Wigner functions are different. Once again, the symmetrical order is 
extraordinarily sensitive to any fingerprint of nonclassical behavior. Without 
going into mathematical details, wc merely quote that, as shown in Ref. [52], 
the SU(2) Wigner kernel can be represented as 

*s°u(2)(^,'/') = U{9,<P) F{1) U\e,4>). (33) 



where F[Jz) is an operator whose expression is of little interest for our pur- 
poses here. It turns out, however, that in the asymptotic case of large dimen- 
sions of the representation (J » 1), F[Jz) tends to the parity operator on 
the sphere: F{Jz) — > exp(— iTrJ^). The SU(2) Wigner function coincides thus 
with (32), except for the factor the appears premultiplying, which does not 
play any relevant role. So, in the classical limit, both approaches give the same 
result. Nevertheless, we mention that the SU(2) quasidistributions have been 
not determined experimentally yet, at difference of the easy measurability of 
the s-ordered quasidistributions for the single polarization modes. 

We finally note that the operator (32), when restricted to a single SU(2) 
invariant subspace, does not contain complete information: in other words, it 
cannot be inverted to obtain the density matrix. 
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4 Examples and concluding remEirks 



To gain further insights into this formahsm, we shall particularize the Wigner 
function (32) for several states of interest. For simplicity, we assume a pure 
state 1^), so that the Wigner function is simply 

W{e,<t^) = Tr[g{e,<j>) {-lf-'n2aUH + i)], (34) 



where 

g{e,ci>)^u\e,ci>) |*)(*| u{e,ci>). (35) 

We recall that, for fixed A^, the states \N, k) {k = 0, . . . , N) span a (A^ + 1)- 
dimensional invariant subspace, wherein the action of the operators (29) is 
standard. In consequence, we expand |^') in this basis 

oo A'' 

EJ2'^Nk\N,k), (36) 

N=0 k=0 

so that the action of U{9,4>) on |^) can be easily calculated. The resulting 
Wigner function takes the form 

oo TV 

Wi9,<f>)=Y: E (-l)"(2iV-2n+l) 

N=0 k,m,n=0 

X *iv.*^^ e^-^^^-) C.(-^)Cn(^), (37) 
where d^^{6) is the Wigner d function 

d^,^{e) ^ {N,k'\eMi&Jy)\N.k), (38) 

whose properties have been extensively studied [53]. Using such properties and 
after some lengthy, but otherwise straightforward calculations, we finally get 

oo 

W{e,<P)=Y, J2 {-IT[N + l + {m + n- N)/cos0] 

N=Q m,n=0 

x*^,„*;,„e'*(™-")dL(-2^)- (39) 
It can be checked that W{9, (j)) is properly normalized 

^ j dnw{e,4>)^i. (40) 
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We first consider the case of (quadrature) coherent states 



l"^) = \aH)H®\av)v- (41) 



We take both oscillators with the same amplitude and relative phase (/?; i.e., 
an = r/ \pl and ay = r/\/2e~'"^. This is perhaps the most interesting situation 
as polarization is concerned. The decomposition in invariant subspaces reads 

= e-l"l^ I I e'^^ (42) 



Note that this state is separable and ^Nk is just the product of the photon- 
number amplitudes in each polarization mode. The sums in (39) can be carried 
out explicitly, with the result 



W{9,(l)) = {r^[l + sinecos(0 + (^)] + l} exp {-2r^[l - sin ecos(0 + (/?)]} . 

(43) 



To examine how the quantum character of the state reflects itself in the prop- 
erties of the Wigner function, we plot the distribution [28] 

(44, 



on the unit sphere. Notice that, with our parametrization, if the unit vector n 
points in the positive direction x, y, or z the state is right-circularly, linearly at 
45°, or horizontally polarized light, respectively. In Fig. 1 we have plotted the 
normalized probability distribution f{6, (f)) for the coherent case, with r — 5 
and the relative phase ip = 7r/2, which classically corresponds to circularly 
polarized light. What we see in the figure is indeed a Gaussian distribution 
centered at the corresponding classical point. 

Next we consider coherent squeezed states in both polarization modes. In 
consequence, we have 

l"^) = \aH,^H)H^\av,^v)v, (45) 



where 

= %)!«)• (46) 
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Fig. 1. Spherical plot of the normalized distribution (44) when the two polarization 
modes are in a coherent state. We take both amplitudes equal to r = 5 and the 
relative phase ip = tt/2. The Gaussian is centered at the classical point corresponding 
to right-circularly polarized light. 

Here \a) denotes a coherent state and 



S{^) = exp 



1 



t2^ 



(47) 



is the squeeze operator. The complex number ^ is known as the squeezing 
parameter and is usually expressed ^ = |.^|e*'^: while |^| measures the squeezing 
of the fluctuations, d measures the direction in which such a squeezing takes 
place. 

The state is again separable and the decomposition (36) is now 

^Nk = h{N- k\aH, ^h)h v{k\av, ^v)v- (48) 



The photon-number amplitude for each polarization mode is 



{k\a,0 = ^ exp 



— a 
2 I' ' 



-a 



(49) 



where B.k{x) are the Hermite polynomials and we have used the notation 
/i = cosh|e|, u = e'^smh\^\. (50) 
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Fig. 2. Spherical plot of the normalized distribution (44) when the two polarization 
modes are in a squeezed state. We take both modes in the same state, with a 
coherent amplitude a = 5 and a squeezing factor |^| = 0.3 in the direction = 0. 

In this case, we have found no simple closed expression for the Wigner func- 
tion (39). Nevertheless, numerical calculations are simple. In Fig. 2 we have 
represented the normalized distribution (44) when both modes are in the same 
squeezed state with a real coherent amplitude a = 5 and a squeezing factor 
I = 0.3 in the direction = 0. Apart from a small deformation of the sphere, 
we see two symmetrical peaks with elliptical contours that indeed represent 
squeezing of the fluctuations. 

To give another simple but illustrative example, we consider a two-mode 
squeezed vacuum state 



and fi and u are defined as in equation (50). The decomposition in invariant 
subspaces for this state is 



^) = smo)H^\o)v), 



(51) 



where the two-mode squeeze operator is 




(52) 




(53) 
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Fig. 3. Spherical plot of the normalized distribution (44) for a two-mode squeezed 
vacuum state with a squeezing parameter |^| = 0.9. 

The Wigner function (39) can be expressed again in a closed form: 



In Fig. 3 we have plotted the normalized distribution corresponding to a two- 
mode squeezed vacuum with |^| = 0.9. We see the presence of two Gaussian 
cups centered at the north and south poles and also a belt around the equa- 
tor of the unit sphere. This state can be seen as arising mainly from these 
contributions and is rotationally symmetric, as (54) is independent of (f). 

Finally, we consider the propagation of light in a Kerr medium. If initially 
both polarization modes are in coherent states of amplitudes an and ay-, the 
state at time t can be written as [54,55,56,57,58] 

\^{t)) =f{T){\aH)H®Wv)v), (55) 



where 

T(r) = exp (ird^fjan aj/flv) , (56) 



and r = x^, X being a real parameter proportional to the third-order nonlinear 
susceptibility of the medium. The state (55) is entangled, and its expansion 
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Fig. 4. Spherical plot of the normalized distribution (44) for light propagating in 
a Kerr medium. Initially the state is coherent in both polarization modes with 
amplitude r = 5 and the (adimensional) propagation time is r = 7r/2. 



(36) can be easily worked out 



The presence of the quadratic phase in (57) prevents again from obtaining an 
analytical form for the Wigner function. In Fig. 4 we show the normalized 
distribution on the sphere for a coherent amplitude r = 5 equal in both po- 
larization modes and r = 71/2. The complicated phase dynamics predicted by 
the theory manifests itself in a complex pattern of well resolved peaks. 

In summary, what we expect to have accomplished in this paper is to present 
a simple alternative phase-space formalism for polarization on the Poincare 
sphere, based on s-ordered quasidistributions for the two basic polarization 
modes. These results may have interesting experimental consequences in order 
to implement a feasible experimental procedure for determining polarization 
properties. 



.iT(N-k)k g-(|aj,|2 + |avf ) 




(57) 
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